Free energy of alternating two- component polymer brushes on cylindrical templates 
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We use computer simulations to investigate the stability of a two-component polymer brush 
de-mixing on a curved template into phases of different morphological properties. It has been 
previously shown via molecular dynamics simulations that immiscible chains having different length 
and anchored to a cylindrical template will phase separate into stripes of different widths oriented 
perpendicularly to the cylindrical axis. We calculate free energy differences for a variety of stripe 
widths, and extract simple relationships between the sizes of the two polymers, A^i and N2, and the 
free energy dependence on the stripe width. We explain these relationships using simple physical 
arguments based upon previous theoretical work on the free energy of polymer brushes. 



Introduction 

Polymer brushes are highly tunable systems that have 
recently attracted considerable attention in the scientific 
community. This is mainly due to the already large num- 
ber of technological applications in which they are used, 
but also to the promising role these materials hold for 
the future. Apart from their well known role in the sta- 
bilization of colloidal particles [Ij, polymer brushes are 
also used as lubricants, in chromatographic devices and 
in adhesives, and their use has recently been proposed in 
a variety of biotechnological applications including drug 
delivery and drug-biocompatibility enhancers [IHS] . Gen- 
erally speaking, they offer an ideal platform that provides 
control over the physical and chemical properties of solid 
and fluid surfaces. 

Polymer brushes are basically dense systems of poly- 
mer chains having one end tethered to a non-adsorbing 
surface, and they have been thoroughly studied theoret- 
ically and numerically on different geometries (see P4T5] 
and references therein). The equilibrium properties of 
these systems typically depend on the molecular weight 
of the chains, their grafting density, and the quality of 
the solvent. Although we have a good qualitative under- 
standing of homogeneous and single components brushes, 
for which scaling arguments have been successfully put 
forward, the problem becomes very complex as soon 
as we deal with non-homogeneous systems with multi- 
component or nonlinear chains. [T6l [TT] 

Of particular relevance for the present paper is the 
work done by Stellacci and collaborators p!8H22] where 
basically a two-component brush of immiscible ligands 
having different length has been shown to phase separate 
into striped phases of different widths when anchored to 
spherical and cylindrical templates. What is surprising 
is not the phase separation in itself, but rather the for- 
mation of the striped phase and its not well understood 
dependence on the mismatch in polymer length. Mi- 
crophase separation is also expected in more traditional 
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mismatched two-component polymer-brushes [23^ , where 
potential applications in stimuli-responsive systems and 
nanotemplating have recently attracted great attention 
to the field (see ref ^ and references therein). Being 
able to control pattern formation at the nanoscale is a 
core element in the production of novel materials via the 
process of self-assembly. 

Direct coarse grained molecular dynamics simulations 
by Glotzer et al. have clearly indicated the driving forces 
behind the formation of these phases as the difference be- 
tween the chains molecular weight, and the overall height 
of the brush itself [H EH EH EH [26] . Very recently, a 
coarse-grained model for two-component polymer mix- 
tures of different length has also been put forward, and 
the formation of different phases as a function of length 
mismatch on a planar geometry has been analyzed, re- 
vealing the presence of striped phases in these systems 
as well [23 ^. 

Unlike other numerical work that focused on the details 
of the de-mixing transition, in this paper we focus ex- 
clusively on the origin of the microphase separation. We 
use numerical simulations to compute free energies of an 
inhomogeneous two-component polymer brush to under- 
stand how molecular weight and overall brush height af- 
fect the stability and width of the striped phases. Specif- 
ically, we focus on the case of a cylindrical solid template 
where stripes have been shown to form promptly both nu- 
merically and experimentally [^ES], and we compute, 
at constant grafting density and template radius, how 
the system free energy changes when imposing stripes 
of different width and height on the template. Our 
results indicate that what limits macrophase separation 
into a two-phase region, thus leading to the stabilization 
of the striped phase, is the elastic strain that builds up in 
the polymer brush formed by the mismatched (exposed) 
polymer segments. 

Methods 

To ensure that our data are not affected by the lateral 
tension due to the immiscibility of the two polymer types, 
and that we are indeed measuring exclusively the free en- 
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FIG. 1: Snapshots of typical initial configurations for 
Lp = 6 and = 30, respectively. 



ergy difference arising from the chains' length mismatch 
and overall brush morphology, we simulate systems in 
which the only difference between the two polymer types 
is their length. This is equivalent to setting the line ten- 
sion between the stripe's boundaries that comes directly 
from the immiscibility term in the pair potential between 
the two polymer types equal to zero. As a consequence, 
we are required to lock in place the position of the an- 
choring monomers to prevent the different chains from 
trivially mixing. 

Polymers are modeled as sequences of spherical beads 
of radius a linearly connected via a harmonic potential 

^ond(^) = ^ ia ~ ^) ' with = 800/cbT. Any two 
monomers in the system interact via the soft and purely 
repulsive dissapative particle dynamics (DPD) simula- 
tion potential: 



LAMMPS [27] with Langevin dynamics in the NVT 
ensemble. Dimensionless units are used throughout this 
paper. The timestep size was set to dt = 0.0025ro (tq is 
the dimensionless time). 

Our goal is to compute the free energy difference 
AF{Lp) between systems having different stripe width 
Lp, while keeping everything else unaltered, as a func- 
tion of the difference between the lengths Ni and N2 . To 
compute the system free energies we use the thermody- 
namic integration method [28 . The idea is to introduce 
a fictitious potential 



if r < <j 
otherwise 



which acts between any two monomers in the system. 
For A = the chains are ideal and the system free energy 
is independent of the specific grouping of the polymer 
types; as A ^ 1 we recover the system of interest. The 
free energy of the full system can then be extracted by 
performing the following integral 



dX 



(1) 



In practice, we perform simulations for several values of 
A, and numerically compute the integral above. 
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where i and j indicate the identity of the polymer i,j G 
{1,2} and m^n e {l^Ni} refers to the identity of the 
monomer. Here Ni is the length of a polymer of type i. 
All the results that follow use en = 622 = ^12 = ^ksT. 
Each polymer is grafted to a fixed point on the outer 
surface of a cylindrical template of radius R = 2.5a 
via the same harmonic potential tethering the consec- 
utive monomers in a chain. In all our simulations we 
considered a total number of polymers Up = 2580 ar- 
ranged in a homogenous grid having square symmetry. 
The polymer identity is finally selected to generate alter- 
nating stripes of width Lp. For a lateral grafting den- 
sity of the polymers equal to 2.74/cr, our system contains 
60 one-polymer-wide rings, making the possible values 
of Lp G {1,2,3,5,6,10,15,30}, when equal number of 



polymer types are considered rip 
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ure [T] shows snapshots of typical initial configurations for 
Lp = 6 and Lp = 30. This particular geometry is 
selected because we find that test unconstrained simula- 
tions of immiscible chains always lead to stripe formation 
perpendicular to the cylindrical axis. This result has also 
been observed in [26 . 

Given the size of the system we carried out our 
simulations using the molecular dynamics package 



Results 

To get insight into the problem we first consider the 
case in which N2 = 0; that is, only one chain type is 
present in the system. Figure |2a| shows our numerical 
data for the free energy cost F{Lp) associated with the 
partitioning of the two polymer species into alternate 
stripes of different Lp on a cylinder of radius 2.5cr and 
lateral grafting density 2.74/ cr^ as a function of Lp and 
for polymer lengths Ni G [10, 30] and N2 = 0. This spe- 
cific lateral grafting density was selected because uncon- 
strained simulations under this condition lead to prompt 
microphase separation for moderate values of polymer 
immiscibility. Furthermore, a sufficiently large value of 
lateral grafting density guarantees large enough differ- 
ences in free energies as a function length mismatch. 

Our data show that the most favorable state corre- 
sponds to that having the largest number of stripes (i.e. 
the smallest value of Lp). The inset on the same figure 
shows how the free energy difference AF^^^ = F{Lp = 
30) — F{Lp = 1) grows with the length of the chain A^i. 
These results can be qualitatively understood by referring 
back to the theory of polymer brushes on fiat surfaces and 
its extension to spherical and cylindrical surfaces O [29] . 
The key point is that the main contribution to the free 
energy cost per chain associated with a polymer brush 
has an inverse dependence on the the lateral chain-to- 
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(a) I Free energy per polymer F/up vs. Lp for a 
inder of radius 2.5<j, for various values of 
Ni e [10, 30] with N2 = 0. Each line represents a 
different value of Ni; increasing Ni leads to an increase 
in free energy F at every value of Lp. The inset shows 

how the free energy gap AF^^^ defined as the free 
energy difference between Lp = 30 and Lp = 1 depends 
on the overall length of the chains Ni. |(b)|The same 



data plotted as a 
chain grafting distance D. 
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where ho ^ crND~^^^ is the height of a brush made with 
polymers of length TV grafted on a planar surface. Al- 
though in our system the location of the grafting lattice 
on the cylindrical surface is independent of Lp and there- 
fore the total grafting density of the system is kept con- 
stant at all times, different distributions of the chains 
on the lattice result in different local densities. Cru- 
cially, when considering the free energy contribution due 
to the different lateral organization of the chains, one 
can consider our system with = 1 as a polymer brush 
with twice the lateral grafting distance of that relative to 
Lp — 30. Clearly the latter system cannot be straightfor- 
wardly described with Eq.[2]due to the free boundaries of 
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FIG. 3: Sketch of expected chain configurations for 
I/„ = 1, 2, and 4. 



the brush, nevertheless we expect that for Lp ^ 00 the 
role of the boundary should become negligible. If we now 
introduce the dimensionless parameter a{Lp) defined as 



FjLp) - F(l) 
F(l) ' 



(3) 



the argument above leads to clear upper-bounds on the 
value of a. Specifically, we expect amax = 2^/^ — 1 in 
the fiat limit, and a^ax = 2^/"^ - 1 for {ho/R) > 1, 
independently of chain length. Our numerical results are 
indeed within these bounds. 

Figure [2b] shows a plotted as a function of Lp for dif- 
ferent values of Ni. a has a nontrivial dependence on 
Lp; nevertheless, simple geometrical arguments can be 
used to understand at least qualitatively the overall be- 
havior of these curves. Figure [3] shows a sketch of the 
expected chain distributions for Lp = 1,2 and 4 for the 
simpler case of a fiat surface. For small values of Lp the 
problem is dominated by boundary effects, and we can 
imagine the chains equally sharing the overall space 2Lp 
per stripe available to them. 

It is immediately obvious that the free energy differ- 
ence between Lp = 1 and Lp = 2 per chain should be very 
small indeed. In both cases the chains can be accommo- 
dated into space as a sequence of blobs of diameter equal 
to twice the grafting distance D. A similar argument can 
be made for Lp = 4. The main extra cost in free energy 
per polymer for a growing Lp should come from the weak 
entropic stretching of the polymers next to the grafting 
surface and should weakly grow as L^. 

When Lp becomes large, the free energy of the system 
is dominated by chains in the bulk of the stripe. In this 
regime the blob size of the outer chains is larger than 
that of the inner chains. Assuming a generic size pro- 



file for the blobs of the form D{i) = D{1 



) where 



i G [l,Lp] and c is fixed by the constraint X^^^ D{i) = 
2DoLp^ one can estimate the free energy per stripe as 
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FIG. 4: a vs. the rescaled Lp/hB- The black hne is a 
guide to the eye. 



FIG. 5: a vs. the rescaled Lp/hs for different values of 
N2. The lines are guides to the eye. 
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Up/Lp particles in each stripe, the free energy per particle 
on a fiat interface reads as 



E 



(4) 



It is difficult to extract P from our simulations, but from 
visual inspection of the brush profiles one should expect 
/3 > 2, leading to functional forms for F having the typ- 
ical saturating behavior shown in our data. 

In the current form the data do not collapse onto a uni- 
versal master curve. In fact, the free energy difference 
of shorter chains tends to grow faster than that of the 
longer ones. This result is quite revealing as it directly 
reflects on the nontrivial interactions between neighbor- 
ing regions of grafted polymers. 

Such interactions will cease as soon as the mutual dis- 
tance between two neighboring regions is larger than 
twice the side spread of the polymers at the boundaries. 
Clearly this will happen sooner for the short polymer 
than for the longer ones. We should be able to account 
for this difference by rescaling of Lp with Lp/N^ . We 
find the best collapse for 7 = 0.4. This result suggests a 
length-scale, hs^ controlling the extent of the side inter- 
actions that scales as /ib = o-ho/Rc' Figure H shows the 
data collapse into a universal master curve a = <l>(^). 

Let's now turn our attention to the case in which 
A^2 7^ 0. It would be desirable to develop a simple way of 
extending our results to the more general N2 ^ case. 
In particular, we are interested in understanding whether 
it is possible to map a system having different polymer 
lengths A^i > A/'2 7^ to a system of a single polymer 
type of length N[ = Ni - N2 (and = 0.) 

Figure [5] shows a as a function of Lp/hB for differ- 
ent values of A^i at A/'2 = 5 and N2 = 10. The results for 
N2 = are also shown as a reference. As the cause for the 
free energy gap between the different configurations rests 
in the difference Ni — A^2, to properly compare the data 



coming from systems having different values of N2 , when 
computing a we subtracted from the measured free ener- 
gies F{Lp) the constant core free energy of the underly- 
ing full brush of height N2, Fcore(A^2). When N2 ^ 0, we 
define a as a{Lp) = {F{Lp) - F(1))/(F(1) - F,,,,{N2)). 

It is comforting to report that indeed it is possible 
to systematically collapse all of the data correspond- 
ing to different values of Ni at a given (fixed) N2 into 
unique curves; however, we cannot collapse data coming 
from different values of A^2- This result clearly points 
to the fact that the universal curve previously described 
for N2 = also has a nontrivial dependence on N2: 
4> = <l>(^, A/'2). In principle, this shouldn't come as a 
surprise ifor brushes grafted on a curved surface, in fact, 
the effect of the "interior" layer of N2 monomers is to in- 
crease the effective radius of the cylinder to which the re- 
duced system would be grafted onto, and to decrease the 
grafting density of the effective brush of length A^i — N2] 
both are variables the free energy depends on (see Eq. [2| . 

Simple geometrical arguments can be used to estimate 
the change in lateral grafting distance. Following [29] . 
we can write the change in lateral grafting spacing as 
we radially move away from the cylinder surface of an 
amount r as ^(r) = D{r / R)^^'^ . For short polymers and 
at the relatively high densities of our systems we have 
r R^ (7N2. It follows that ^(r) = D{1 + aN2/Ry^'^. 
Simultaneously, the effective polymer of length A^i — 7V2 
would be grafted onto a cylinder of radius R^ R-\- (JN2 . 
The ratios R' appearing in the brush free energy should 
then be adjusted to R' /i{r) {R/D){1 + aN2/Ry^^ 
introducing a predictable dependence on N2 in the prob- 
lem. Although this argument suggests an overall increase 
in free energy for the system at a given Lp, it isn't obvious 
how it will affect the dimensionless ratio a{Lp). Clearly 
this correction should factor out in the large R limit or 
for ho/R > 1. 

Our numerical data show that a is quite sensitive to 
N2. Specifically, a{Lp) tends to flatten as N2 increases. 
This is most likely explained by the fact that the base of 
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the effective chain of length Ni — N2 is not reahy grafted 
in place, but the layer of polymers below allow for signifi- 
cant reorganization of the base due, for instance, to their 
compressibility. Because of the blob-size dependence on 
the distance from the surface, on a convex template the 
chains become systematically more laterally compressible 
as N2 is increased. The net result is a flattening of the 
free energy difference as a function of Lp. Unfortunately, 
we find that our data for different values of A^2 cannot be 
collapsed by a simply rescaling of the free energy with a 
power of indicating a subtle interplay between F{Lp) 
and the height of the bottom layer. 

Conclusions 

In summary, we analyzed the relative stability of the 
striped phases arising from the de-mixing of immiscible 
polymers on a cylindrical template as a function of the 
width of the stripe, the length mismatch between the 
chains, and the overall length of the brush. We set nu- 
merical bounds for the free energy gap between an al- 
ternating (many thin stripes) and a fully de-mixed (two 
wide stripes) phase, and found that indeed the former be- 
comes more and more favorable as the mismatch between 
the chains' lengths, AA^ is increased. We also found 
that when appropriately normalized the free energies as a 
function of stripe width for different values of can be 
collapsed into unique master curves that only depend on 



the overall length of the lower brush N2. Finally we dis- 
cussed the possibility of reducing a two-component sys- 
tem of chains having lengths Ni > N2 into a simpler 
single component system of chains of length AN. Our 
results are fully consistent with previous molecular dy- 
namics simulations on this problem in identifying the 
key parameters setting the width of the stripes in the 
polymer length mismatch N1—N2, and the overall length 
of the shorter brush 7V2 together with the degree of chain 
immiscibility. 

In all our data the direct line tension 7 coming from the 
pair potential between chains of different types has been 
set to zero. This term would add to the free energy bal- 
ance a predictable contribution of the type Fi{Lp, N2) ^ 
7 hQ{N2){Np/Lp)^ where ho{N2) is the height of the poly- 
mer brush formed by the lower chains and {Np/Lp) is 
the number of phase boundaries. Clearly Fi{Lp^N2) has 
a minimum for Lp ^ 00 and is large for small Lp. The 
balance between the line tension and the configurational 
free energy computed in this paper should therefore set 
the width of the stripe. More work in this direction is 
currently underway. 
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